CENTRAL LIMIT THEOREMS FOR THE BROWNIAN MOTION ON 

LARGE UNITARY GROUPS 



FLORENT BENAYCH-GEORGES 



Abstract. In this paper, we are concerned with the large n hmit of hnear combinations 
of the entries of a Brownian motion on the group of n x n unitary matrices. We prove that 
the process of such a hnear combination converges to a Gaussian one. Various scales of 
time and initial distribution are concerned, giving rise to various limit processes, related 
to the geometric construction of the unitary Brownian motion. As an application, we 
propose a quite short proof of the asymptotic Gaussian feature of the linear combinations 
of the entries of Haar distributed random unitary matrices, a result already proved by 
Diaconis et al. 



Introduction 

There is a natural definition of Brownian motion on any compact Lie group, wliose 
distribution is sometimes called the heat kernel measure. Mainly due to its relations with 
the object from free probability theory called the free unitary Brownian motion and with 
the two-dimentional Yang- Mills theory, the Brownian motion on large unitary groups has 
appeared in several papers during the last decade. Rains, in |R97] . Xu, in |X97] . Biane, in 
|B97a| IB97b] and Levy and Mai'da, in |L08| ILM09] , are all concerned with the asymptotics 
of the spectral distribution of large random matrices distributed according to the heat 
kernel measure. Also, in |D08j . Demni makes use of the unitary Brownian motion in 
the study of Jacobi processes, and, in |BL08j . Levy and the author of the present paper 
construct a continuum of convolutions between the classical and free ones based on the 
conjugation of random matrices with a unitary Brownian motion. In this paper, we are 
concerned with the asymptotic distributions of linear combinations of the entries of an 
n X n unitary Brownian motion as n tends to infinity. 

We first give the joint limit distribution, as n tends to infinity, of (possibly rescaled) 
random processes of the type {TT[A{Vt — I)])t>o for (e~*''^VJ)t>o a Brownian motion staring 
at / on the group of unitary n x n matrices and A a random n x n matrix, independent 
of {yt)t>o (Theorem 11.21) . This theorem is the key result of the paper, since specifying 
the choice of the matrices A will then allow us to prove all other results. As a first 
example, it allows us to find out, for any sequence (a„)„ of positive integers with a limit 
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a G [0,+oo], the limit distribution of any upper-left corner of \/n/a^{V\og{ant+i) — I)t>o 
(Corollary II. 3p : for small scales of time (i.e. when a = 0), the limit process is purely 
skew-Hermitian, whereas for large scales of time {a = +00), the limit process is a standard 
complex matricial Brownian motion. For intermediate scales of time (0 < a < +cxd), the 
limit process is an interpolation between these extreme cases. The existence of these three 
asymptotic regimes can be explained by the fact that the unitary Brownian motion is the 
"wrapping", on the unitary group, of a Brownian motion on the tangent space of this 
group at I (which is the space of skew-Hermitian matrices), and that as the time goes to 
infinity, its distribution tends to the Haar measure (for which the upper-left corners are 
asymptotically distributed as standard complex Gaussian random matrices). 

Secondly, we consider a unitary Brownian motion {e~^^^Vt)t>o which initial distribution 
is the uniform measure on the group of permutation matrices: its entries are exchangeable 
at any value of t. In this case, for any positive sequence (a^) and any positive integer 
p, the p X p upper left corner of (A/n/a„V^og(a„t+i))t>o converges to a standard complex 
matricial Brownian motion (Corollary ll.7p . 

Since the unitary Brownian motion distributed according to the Haar measure at time 
zero has a stationary distribution, our results allow us to give very short proofs of some 
well-known results of Diaconis et al, first proved in |ADN03j , about the asymptotic nor- 
mality of linear combinations of the entries of uniform random unitary matrices (Theorem 
Ol and Corollary fLTOj) . 

Let us now present briefly what problems underlie the asymptotics of linear combina- 
tions of the entries of a unitary Brownian motion. 

Asymptotic normality of random unit vectors and unitary matrices: The historical first 
result in this direction is due to Emile Borel, who proved a century ago, in [B1906] . 
that, for a uniformly distributed point (Xi, . . . ,X„) on the unit euclidian sphere S"'"^, 
the scaled first coordinate y/nXi converges weakly to the standard Gaussian distribution 
as the dimension n tends to infinity. As explained in the introduction of the paper 
|ADNn3] of Diaconis et al, this says that the features of the "microcanonical" ensemble 
in a certain model for statistical mecanics (uniform measure on the sphere) are captured by 
the "canonical" ensemble (Gaussian measure). Since then, a long list of further- reaching 
results about the entries of uniformly distributed random orthogonal or unitary matrices 
have been obtained. The most recent ones are the previously cited paper of Diaconis et 
al, the papers of Meckes and Chatterjee [MO 8 1 ICMOSj . the paper of Collins and Stolz 
|CS08j and the paper of Jiang |J06j . where the point of view is slightly different. In the 
present paper, we give a new, quite short, proof of the asymptotic normality of the linear 
combinations of the entries of uniformly distributed random unitary matrices, but we also 
extend these investigations to the case where the distribution of the matrices is not the 
Haar measure but the heat kernel measure, with any initial distribution and any rescaling 
of the time. 

Second order freeness: A theory has been developed these last five years about Gaussian 
fluctuations (called second order limits) of traces of large random matrices around their 
limits, the most emblematic articles in this theory being |MNn4[ iMSOHl IMSSn7[ KMSSO?] . 
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The results of this paper can be related to this theory, even though, technically speaking, 
we do not consider the powers of the matrices her^. 

Brownian motion on the Lie algebra and ltd map: The unitary Brownian motion is a 
continuous random process taking values on the unitary group, which has independent 
and stationary multiplicative increments. The most constructive way to define it is to 
consider a standard Brownian motion {Bt)t>o on the tangent space of the unitary group at 
the identity matrix and to consider its image by the ltd map (whose inverse is sometimes 
called the Cartan application), i.e. to wrap@ it around the unitary group: the process 
{Ut)t>o obtained is a unitary Brownian motion starting at I. Our results give us an idea 
of the way the Ito map adulterates the Brownian feature of the entries of {Bt)t>o at 
different scales of time. Moreover, the question of the choice of a rescaling of the time 
(depending on the dimension) raises interesting questions (see Remark 11.11) . 

Notation. For each n > 1, U„ shall denote the group of n x n unitary matrices. The 
identity matrix will always be denoted by /. For each complex matrix M, M* will denote 
the adjoint of M. We shall call a standard complex Brownian motion a complex-valued 
process whose real and imaginary parts are independent standard real Brownian motions 
divided by \/2. For all > 1, the space of continuous functions from [0, +oo) to C*^ 
will be denoted by C([0, +oo), C'^) and will be endowed with the topology of the uniform 
convergence on every compact interval. 



1.1. Brief presentation of the Brownian motion on the unitary group. There 
are several ways to construct the Brownian motion on the unitary groupy. For the one 
we choose here, all facts can easily be recovered by the use of the matricial Ito calculus, 
as exposed in Section 12.11 

Let n be a positive integer and z/q a probability measure on the group of unitary n x n 
matrices. We shall call a unitary Brownian motion with initial law uq any random process 
{Ut)t>o with values on the space of n x n complex matrices such that Uq is i^o-distributed 
and {Ut)t>o is a strong solution of the stochastic differential equation 



The reason is that the constant matrices we consider here, hke ^/nx (an elementary n x n matrix), 
have no bounded moments of order higher than two: our results are the best ones that one could obtain 
with such matrices. 

^For G a matricial Lie group with tangent space g at /, the "wrapping" Wj, on G, of a continuous 
and piecewise smooth path 7 : [0,+oo) g such that 7(0) = is defined by Wj{0) — I and w'^{t) = 
w^{t) X 7'(i). It -B is a Brownian motion on g and {Bn)n>i is a sequence of continuous, piecewise afhne 
interpolations of B with a step tending to zero as n tends to infinity, then the sequence wb^ converges in 
probability to a process which doesn't depend on the choice of the interpolations and which is a Brownian 



motion on G pml Sect. VI.7], |RW87l Eq. (35.6)], [FO09] . 

^See [H56l ISV73[ llWSli iRWg?] . A very concise and elementary definition is also given in [R97j . 



1. Statement of the results 
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where {Ht)t>o is a Brownian motioiiB on the space of n x n Hermitian matrices endowed 
with the scalar product {A,B) = nTT{AB). 

It can be proved that for such a process {Ut)t>o, for any ^ 0, f/^g is almost surely 
unitary and both processes {UtQ+tU^^)t>o and {U^^Uta+t)t>o are unitary Brownian motions 
starting at /„ and independent of the cr-algebra generated by {Us)o<s<to- 

Remark 1.1 (Communicated by Thierry Levy). There are other ways to scale the time 
for the Brownian motion on the unitary group. Our scaling of the time is the one for 
which the three limit regimes correspond respectively to small values of t, finite values of 
t and large values of t and for which the limit non commutative distribution of {Ut)t>o is 
the one of a free unitary Brownian motion. It also has a heuristic geometrical meaning: 
with this scaling, for any fixed t, the distanc^ between Uq and Ut has the same order as 
the diameter^ of the group. It means that for any fixed t > 0, large values of n, Ut is 
probably no longer too close to its departure point, while it also probably hasn't "orbited" 
the unitary group too many times. 

1.2. The three asymptotic regimes for the unitary Brownian motion starting 

at In- Let (an)n>i be a sequence of positive numbers with a limit a G [0,+oo]. Let us 
fix a positive integer k and let, for each n > 1, A", . . . , he a family of random n x n 
matrices and (e~*/^VJ")t>o be a Brownian motion on U„ starting at J, independent of 
the y4"'s. Suppose that there exists nonrandom complex matrices [ai]f=i, \pi,i']fi'=i and 
[Qi,i']f,i'=i such that we have 



(2) 


V/, 


-Tr(Ar) ai 

n n— s-oo 


(converj 


^ence 


in probability). 


(3) 


v/,r. 


-Tt{A^A^,) 

n n— >-oo 


(converj 


^ence 


in probability). 


(4) 


v/,r. 




(converj 


^ence 


in probability). 


For 


each n, let 


US define, for t > 0, 









= a-'/'iTT[A^M:^^,^t+i) -/)],•••, nAi{Vy:,^a„t+i) -m- 

Let /i be the probability measure on C([0, +oo), C*^) defined by the fact that any process 
Xt = {Xt^i, . . . , Xt^k) distributed according to yU is a Gaussian centered process with 



{iHt)t>o is in fact the skew-Hermitian Brownian motion that the process (?7t[/o)t>o wraps around 
the unitary group, as explained in Footnote [51 

^As explained in Footnotes [5] and 01 (Ut)t>o is the wrapping of {iHt)t>o on the unitary group, hence 
the distance between Uq and Ut has the same order as the one between and iHt, which, by the Law of 
Large Numbers, has order \/tn. 

^By definition, the diameter of the group is the supremum, over pairs U, V of unitary matrices, of 
the length of the shortest geodesic between U and V. Here, it is equal the length of the geodesic 
t £ [0, n7r] i~> exp{itln/n) between /„ and — /„, i.e. to rnr. 
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independent increments such that for all t > 0, for all 1,1' = 1, . . . , k, 

-PLi't if a = 0, 

(5) E(X,zX^) = q^d , E(X,zX,,,) = { -pi/-^^^^ + a^a/^^^^ if < a < +oo, 

if a = +00. 

Theorem 1.2. As n tends to infinity, the distribution of the process (X")t>o converges 
weakly to fi. 

Recall that a principal suhmatrix of a matrix is a matrix obtained by removing some 
columns, and the rows with the same indices. 

Corollary 1.3. Let us fix p > 1 and let (Ht), {St) be two independent standard Brown- 
ian motions on the euclidian spaces of p x p respectively Hermitian and skew-Hermitian 
matrices endowed with the respective scalar products {X,Y) = Tt{XY)/2, {X,Y) = 
— Tr{XY)/2. Then, as n tends to infinity, the distribution of the VP^"^ -valued process 
of the entries of any p x p principal submatrix of \AV^(^og(a„t+i) ~ ^)t>o converges to 
the one of the random process {Ht-^f^(t) + St+f^{t))t>o, where 

t if a = 0, 

fa{t) = {'-^^^ zfO<a<+oo, 
if OL = +00. 

Remark 1.4. Note that when a = 0, the limit process is simply a standard Brownian 
motion on the space oi p x p skew-Hermitian matrices, whereas, as a grows from zero 
to +00, the Hermitian part of the limit process keeps growing, and at last, when a = 
+00, the limit process is a standard Brownian motion on the space of p x p complex 
matrices. As said in the introduction, the existence of these three asymptotic regimes 
can be explained by the fact that the unitary Brownian motion is the "wrapping", on 
the unitary group, of a Brownian motion on the tangent space at I (which is the space 
of skew-Hermitian matrices), and that as the time goes to infinity, its distribution tends 
to the Haar measure (for which, as stated by Corollary I1.10[ the upper-left corners are 
asymptotically distributed as standard complex Gaussian random matrices). 

1.3. The particular case of unitary Brownian motions with exchangeable en- 
tries. Let (a;„)„>i be a sequence of positive numbers (no hypothesis is made on its con- 
vergence). Let us fix a positive integer k and let, for each n > 1, A^, . . . , he a. family 
of random nx n matrices and (e~*/^V^")j>o be a Brownian motion on U„, independent of 
the A^'s, such that Vq is uniformly distributed on the group of matrices of permutations 
of {1,...,?T,}. Suppose that there exists a nonrandom complex matrix [<ii,v]ui=i a 
probability measure /io on C'^ such that 

(6) V/,/', -Tr[AfAp*] — ^ gn' (conv. in probability), 

(7) V/, -Kiiihi) ; ((^, j)-th entry of A") ^0} — ^ (conv. in probability). 



(8) a-^'^{Ti[AlV^l Tt[AIV^]) — > fiQ (conv. in distribution). 



6 



FLORENT BENAYCH-GEORGES 



For each n, let us define, for t > 0, 

= "n^''^(Tr(A"V^"g(^^j^i)), . . . , Tr(y4^\/i"g(„^4_^i))). 

Let P be an Hermitian matrix sucli tliat = [qi,i']iii=i, {Zi, . . . , Zk) be an independent 
family of standard complex Brownian motions and C be a yUo-distributed random variable, 
independent of the Z;'s. Let us define to be the distribution, on C([0, +00), C'^), of the 
process (C + (Z^^, . . . , Zk,t)P)t>o- 

Theorem 1.5. As n tends to infinity, the distribution of Xn converges weakly to fi. 

Examples 1.6. Let us give a few examples of sequences which satisfy the hypotheses 
©, (IZD and &. 

a) Firstly, if (q;„) is bounded from below by a positive constant, if holds and if 
for all /, n~^^{{i,j) ; ((z, j)-th entry of A") 7^ 0} converges in probability to zero, then ([H]) 
holds for the Dirac mass at zero (this can easily be deduced from Lemma [2.31 bellow). 

b) Secondly, if, a„ tends to one and if, for all n, the matrices A^, 1 < I < k, are random 
real n x n matrices which satisfy and ([7j), then (jS]) holds for /xq the Gaussian measure 
with covariance matrix [qi,v]'ii'=i- This follows easily from |S88] . As an example, if /c = 1 
and if, for each n, the entries of are i.i.d. with distribution r„ such that 

rnm) ^ 1 and n(l - r„({0})) / tMr„(t) 1, 

then ([6]), ([7]) and ([8]) hold for qi^i = 1 and /xq the standard Gaussian law. 

c) Other examples can be found using |078l Th. 5.1], where the laws /iq are other 
infinitely divisible laws. 



Both following corollaries are direct applications of the previous theorem, the first one 
using implicitly the fact that any entry of V^^ is null with probability 1 — n~^, and the 
second one using implicitly the fact that the distribution of the number of fixed points of 
a uniform random permutation of {1, . . . ,n} converges weakly, as n tends to infinity, to 
the Poisson distribution with mean one |DS94 j. 

Corollary 1.7. Let (an) be a sequence of positive integers. For any p >1, as n tends to 
infinity, the distribution of any p x p principal submatrix of 

converges weakly to the one of an independent family of standard Brownian motions on 
the complex plane. 

Corollary 1.8. As n tends to infinity, the distribution of (Tr(V^"gj.j^-|^-j))t>o converges 
weakly to the one of [C + Zt)t>o, where C is a Poisson random variable with mean one 
and {Zt)t>o is a standard complex Brownian motion, independent of C. 
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1.4. Application to the asymptotics of the uniform measure on the unitary 
group. Since the Brownian motion on the unitary group distributed according to the 
Haar measure at time zero has a stationary distribution, our results allow us to recover 
certain results of asymptotic normality of linear combinations of the entries of uniform 
random unitary matrices. 

The following theorem is not really new |CS08t ICM08] , except for the minor improve- 
ment that we allow to "randomize" the matrices A^^s. However, our method allows to 
give a very direct proof, even under these very general hypotheses. 

Theorem 1.9. Let us fix k > 1 and let, for each n > 1, A", . . . ,y4^ be n x n random 
matrices and Un be a random matrix with uniform distribution on the group of n x n 
unitary matrices, independent of the 's. Suppose that for all I, I' , there is a nonrandom 
complex number qi^i such that 

— Tt^A^A^*) — y qui (conv. in distribution). 

n n— >-oo ' 

Then as n tends to infinity, the distribution of the random vector (Tr[y4"t/„], . . . , Tr[A^f/„]) 
converges weakly to the one of a Gaussian centered family {Zi, . . . , Z^) of complex random 
variables such that for all 1,1' = 1, . . . , k, K{ZiZi/) = and ^[ZiZir) = qi^y. 

The following corollary is immediate. 

Corollary 1.10. Let us fix k > 1 and let, for each n, Z^, . . . , ZJ^ be k different entries 
of a random nxn matrix with uniform distribution on the unitary group. Then the joint 
distribution of y/n{Z'^, . . . , Z]^) converges weakly, as n tends to infinity, to the one of a 
family of independent standard complex Gaussian random variables. 

2. Proofs 

2.1. Preliminaries on matricial Ito calculus, a) Let n he a. positive integer. Let 
{J^t)t>o be a filtration and {H^)t>Q be an (J^t)t>o-standard Brownian motion on the space 
oi n X n Hermitian matrices endowed with the scalar product (A, B) = nTr{AB), i.e. a 
process with values in the space of n x n Hermitian matrices such that the diagonal and 
upper diagonal entries of (y^/7")t>o are independent random processes, the ones on the 
diagonal being standard real Brownian motions and the ones above the diagonal being 
standard complex Brownian motions. If one considers two matrix-valued semimartingales 
X, Y such that 

dXt = A{dH^)Bt + Ctdt, dYt = DtidH^Et + F^dt, 
for some (J-i)t>o-adapted matrix-valued processes A, B, C, D, E, F, then, by Ito's formula, 

d{XY)t = {dXt)Yt + XtdYt + -Ti{BtDt)AtEtdt, 

n 

d(Tr(X),Tr(y))i = -Ti{BtAEtDt)dt. 

n 

We shall use these formulas many times in the paper, without citing them every time. 
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b) With the same notation, let us now consider a deterministic function / with 
positive derivative such that /(O) = 0. Then for X the process introduced above, the 
process Xt := -^/(t) satisfies 

dx, = v^A/w(di^r)i?/w + f'it)Cmdt, 

where is the (J-/(t))t>o-Brownian motion defined by the formula HJ^ = J-^ — d(i?/(.))g. 

2.2. Proof of Theorem 11.21 Let us first state some matricial inequalities we shall often 
refer to in the following. Let X, Y be two complex matrices and G, H be two Hermitian 
nonnegative matrices. Then we have 

(9) |Tr(XF)| < v^Tr(XX*)v/Tr(rF*), 

(10) Tr(G2) < (TrG)2, 

(11) |Tr(GiJ)| < Tt{G)Tt{H). 

Inequality ( TTT]) follows from and (ITU]) , which are obvious. 

Lemma 2.1. Let us fix n > 3, an nxn matrix A and a Brownian motion (e~*/^\4)i>o on 
Vn starting at I. Then there exists some real constants Ci, C2, C3, C4 which are bounded 
from above by 100(Tr y4y4*)^/n^ and such that for all t > 0, 

(12) E[Tr{AVtAVt)] = Tt{A^) cosh(t/n) - (Tr A)^ sinh(t/n), 

(13) E[| Tr(AW.)r] = I Tr(A^)P + + {-sh (f ) - l} 

- |3ft[(TrA)2Tr(AA)] + n^/^Ca} sinh 
+nC3(e* - 1) + Ciie^' - 1). 

Proof. Since the formulas we have to state are invariant under multiplication of A by a 
scalar, we can suppose that Tt{AA*) = n. 

Note that for HJ^ as in Section I^TTl (VJ) is a strong solution of dVt = i{dH^ )Vt. Hence 
by the matricial Ito calculus, 

d TiiAVtAVt) = 2i TiiVtAVtAdH]^) - n'^ 1i{AVt) Tr(A\/i)dt, 
d TiiAVt) Ti{AVt) = 2i Tr{AVt) Tr(\4v4di7f ) - n'^ Tr{AVtAVt)dt. 

It follows that for x{t) = E[TT{AVtAVt)] and y{t) = E[TT{AVt) Tr(A\4)], we have 

x' = —n~^y and y' = —n~^x. 

Equation (fT2|) follows. 

Now, let us define, for C,D some nxn matrices, uc,D(t) = lE[Tr(l^CVJ*D)] and 
vcoit) = IE[Tr(VtC) Tr(V^*D)]. By the matricial Ito calculus again, one has 

(14) ucMt) = -(e*-l)TrCTrZ} + Tr(CD), vcMi) = -(e*-l) Tr(CD) + Tr(C) Tr(D). 
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Let US now prove (|T3l) . We introduce the functions 
fit) = E{\TTiAVtAVt)\'), 
g{t) = ^{E[Tr{AVt)TiiAVt)TiiA*V;A*V;)]}, 
h{t) = E[|Tr(Ay,)l1. 

By the matricial Ito calculus again (using the hypothesis Ti{AA*) = n), 

(15) nxf{t) = -2g{t) + Ae'uA*A,AA4t), 
nxg\t) = -f{t)-h{t)+4e'^VA,A*AA*{t)], 
nxh'it) = -2g{t) + Ane%A,A'i^)- 

It follows, by dH, that g"{t) - ^g{t) = ^k + ^9, for 
K = n-^TT{AA*AA*) -1, 

6 = 2 - n-^ Tt{AA* AA*) -n-^Ti{AAA* A*) -iTiAl"^ + ^{TtATt{A*AA*)}, 
hence g{t) = /icosh (f ) +usmh (f ) + ^e^* + ^e*, with 

/. = ^TrATTATT{A*A*)]-^^-4^^ 
u = -^|TrA2|2 - l|TrA|^ + 2gfJ[Trv4Tr(A*AA*)] 
From (1151) . it follows that 

/(t) = -|TrAY-/Usinh(^^^ _ l^cosh (^^^ -1^ +w;(t), 

where w{t) = -j^{e^' - 1) - J^{e' - 1) +2{e'^' - 1) + A{n'^ Tt{AA*AA*) - l)(e*- 

Now, the conclusion follows from the fact that since Tt{AA*) = n, the inequalities ([ 
ffTU]) and (dU) allow to prove that | Ti A\, \ Tr A^j < n, Tt{AA*AA*), \ Tt{AAA*A*)\ < 
and \ Tt{A*AA*)\ < n^/^. 

Lemma 2.2. Let (an), (bn), (c„) be sequences of real numbers such that (an) tends to + 
and (bn) and (c„) are both bounded. Then we have 

(16) := ^ (cosh (i^l - 1 V — sinh fi^") 0. 



— 1 — 4 



Proof. Let us define 

f cosh(x) — 1 sinh(x) 
A = max < sup |a„| , sup |o„| , sup , sup 

[n>l n>l a<x<l X 0<x<l X 

Then f|T6|) follows from the fact that, since Un can also be written 

n%ri n^bn n^br, nCn nCr, 



^" ~ 1-- i+- 1-- 

2an " 2an " 2a„ " 2a 



we have 

I , oi^2log^«n + loga^ 
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□ 

Proof of Theorem 11.21 Suppose the Theorem not to hold. Then there is an infinite 
set M of positive integers, a continuous bonded function / on C([0, +00), C'^) and e > 
such that for all n G M , 

(17) |E(/(X")) - j fM > e. 

Now, recall that if a sequence of real random variables converges in probability, then one 
of its subsequences converges almost surely. Hence one can suppose that the convergences 
([2]), ([3]) and (jl]) hold almost surely as n increases to infinity within the set M . 

We shall prove that as n goes to infinity within the set A/", 



(18) E(/(X")) ^ j /d/x, 

contradicting (ITTI) and closing the proof. 

Up to a conditioning by the a-algebra generated by the y4"'s (1 < / < A;, ?t- > 1), one 
can suppose these matrices to be non-random. 

In this case, for all n, X" is a C^-valued continuous centered martingale. To prove that 
its distribution tends to /x, by ReboUedo's Theorem (see |R80j or |AGZ09[ Th. H.14]), 
it suffices to prove that the bracket of X" converges pointwise, in L^, to the one of 
a //-distributed process. Hence it suffices to fix Ai,...,Afc G C, to define the process 
Yj" = XiX^^ -!-••• + AnXj\ and to prove that as n tends to infinity within the set A/", 




if a = 0, 

(19) (y",F^), ^ qt and (^",1^")^ \ -pM^ + a^ '°g'f if < a < +00, 



2a 



if a = +00. 



for q = Ya,i'=i ^i^i'Qi,i'^ P = Ya,i'=i ^ih'Pi,v and a = XlfLi Aja;. 
Let us define, for each n, A"- = X]f=i '^i^?- We have 

>^" = «n"'/'Tr[A"(V;„V„m)-^)]- 



Hence by Section [2711 b). satisfies 

where is an Hermitian Brownian motion as introduced in Section 12.11 Thus, since 
^og(a„t+i)^og(a„i+i) ~ (^"■^ + ^)^' ^ly the matricial Ito calculus, we have 

d(F",F^)t = -Tr(A"A"*)dl 
n 

Since as n goes to infinity within the set Af, ^ Tt^A^A"^*) tends to q, the first part of ( 1T9|1 
is proved. Let us now prove the second part. By the matricial Ito calculus again, we have 
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Hence it suffice to prove that as n increases to infinity within the set M , we have the 
convergence 

f20) TrW-y" A-V- 1^1^ + "'^^ if < a < +00, 

a = +00, 

uniformly as t varies in any compact subset of [0, +00). This follows easily from Lemma 
12.11 (and Lemma [2.21 in the case where a = +00). □ 

2.3. Proof of corollary II. 3L By Theorem IL2I and formula ([5]), applied with matrices 
Af of the type ^/nx (an elementary n x n matrix), the distribution of any p x p principal 
submatrix of \AV^(^og(a„t+i) ~^)t>o converges weakly to the one of the random process 
(Mi)i>o with independent increments such that for all p x p complex matrices X, Y and 
all t > 0, 

log(at + 1) 



E[Tr(MjX)Tr(MtF)] =tTi{XY*), E[Ti{MtX)Ti{MtY)] = — Tr(Xr), 

a 

where if a = or +00, the second right-hand term has to be replaced by respectively 
— tTr(Xy) or 0. Since a standard Brownian motion (Bf) on an euclidian space {E, (■, ■)) 
satisfies, for all u,v E E, K{{Bt,u){Bt,v)) = {u,v)t, the result can easily be verified. □ 

2.4. Proof of Theorem [Q 

Lemma 2.3. Let us fix n > 2, let S be the matrix of a uniform random permutation of 
{1, . . . ,n} and let A,B he n x n non-random matrices. Then we have 



E{|Tr(A5)|} < -v^/lY(AA 

n 



m\i:iiASBS)\} < - — ^ + '^^f^'' ^ti(aa*)Jtt(bb*), 

n[n — 1) 

where for each matrix X = [xij]^j^i, Cx denotes tl{(^, j) ; Xij 7^ 0}. 

Proof. Let us denote by respectively [aij], [bij], [sij] the entries of A, B, S. By Holder's 
inequality, we have 

E{| Tt{AS)\} < 5^ \a,,M\j) = -Yl I^M'I ^ -V^VTr(AA*). 



n ^ — ' n 



Moreover, we have 



E{\Tt{ASBS)\} < 5^|afc,,-6fe,,|E(4,)+ ^ \a,jbkMsj,kSi,i) 

j,k 

jj^l or k^i 

— ~'y^A^k,jbk,j\ H TT \(li,jbk,l\y 

n ^ nin — 1) ^ 

j,k id,k,l 

and the conclusion follows from Holder's inequality again. □ 

Proof of Theorem 11.51 We consider C, Zi, . . . , Z^, P as introduced above the state- 
ment of the theorem. For each t > 0, let us define Kt = {Zi^t, ■ ■ ■ , Zk,t)P- It suffices to 
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prove that as n tends to infinity, 

(P) the joint distribution of (Xq , (X"— Xg )i>o) converges weakly to the one of (C, {Kt)t>o)- 
Step I. Let us prove the following convergences in probability: 



En:=[yi = l,...,k, Tr(A„,iA;j < n(g,,i + 1) 

and ; ((i,j)-th entry of A^^i) ^ 0} < en^ ] . 



Using Lemma [231 and conditioning the expectation with respect to the cr-algebra generated 
by the An/s, one can easily prove that for all 1,1' = 1, . . . , k, 



Since, by ([6]) and (^^, the probability of the event En tends to one as n tends to infinity, 
the convergences in probability of (12T|) and ( 122|) follow. 

Step II. The weak topology on the set of probability measures on the Polish space 
C'^xC([0, +oo), C'^) is metrizable, thus to prove (P), we have to prove that any subsequence 
of (Xq , (X" — Xq )t>o) admits a subsequence which distribution converges weakly to fi. 
Hence we introduce an infinite set of positive integers A/", and we shall prove that as n goes 
to infinity within the set TV, the distribution of (Xq , (X" — Xq )t>o) admits a subsequence 
converging weakly to the distribution of (C, (i^f)t>o). 

Step III. By Step I and since [0, +oo] is compact, we can suppose that Af is such that as 
n goes to infinity within the set J\f, (an) admits a limit a G [0, +oo] and the convergences 
of (E]), dUD and ([22]) hold almost surely. Note that for all n, (e-*/2-l/Q"*V;")j>o is a unitary 
Brownian motion starting at J, independent of the cr-algebra S„ := cr(A", . . . , A^, Vq). 
By Theorem II. 2[ it implies that almost surely, as n tends to infinity within the set TV, 
the conditional distribution of X" — Xq with respect to S„ converges weakly to the 
distribution of {Kt)t>o- By Hypothesis ([HD, it follows that (P) holds as n increases to 
infinity within A^. It closes the proof of the theorem. □ 

2.5. Proof of Theorem [TS 



Lemma 2.4. Let U be Haar- distributed on and A be an n x n matrix, with n > 3. 
Then 





(23) 
(24) 



E{|Tr(Af/)|2} = TT{AA*)/n, 
E{\Tt{AUAU)\^} < 100(Tr(AA*))V^^. 
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Proof. Set U = [uij]^j^i. One can write A = VA'W, with V,W unitary matrices 

and A' a diagonal matrix whose diagonal entries the eigenvalues of VAA*. 

Since the law of U is invariant under the left and right actions of the unitary group, one 
can suppose that A = A'. These invariances of the law of U also imply that for all i,j, 
K{ui^iUjj) = 51/ n. Equation ( 12^ follows. Equation follows from (IT^ and the fact 
that the Haar measure on the unitary group is the limit of the distribution of e~*/^\4 as 
t tends to infinity. □ 

Proof of Theorem 11.91 Step I. Firstly, following the proofs of Theorems 11.21 and II. 5^ 

we suppose that the y4"'s are nonrandom. Then by Lemma [2. 4^ for all /, Tr(y4"t/„) 
and Ti{A^UnA^,Un) converge in probability to zero as n tends to infinity. 

Step II. Following |DZ98t Th. D8], we fix a bounded real function / on C'^ which is 
1-Lipschitz for the canonical hermitian norm || ■ || on C'^, and we shall prove that 

(25) E{/(Tr[A^f/„], . . . , TrKt/„])} E{/(Zi, . . . , Z,)}. 
Let us fix e > 0. 

a) Let, for each t > 0, {Zi^t, . . . , Z^^t) be a Gaussian family of centered complex random 
variables such that for all I, I', K[Zi^tZi\t] = and K[Zi^tZi'^t] = ~ (such a family 
exists because the matrix [qi^i']i^i' is nonnegative). The distribution of {Zi^t, . . . , Z^^t) tends 
to the one of (Zi, . . . , Z^) as t tends to infinity. Hence there is to > such that 

(26) |E{/(Zi,...,Zfc)}-E{/(Zi,i„...,Zfc,iJ}| < e, 

(27) e-^sup{n-^TrKAr + --- + A^Ar]}'/' < e. 

n>l 

b) For each n, up to an extension of the probability space where [/„ is defined, one can 
suppose that there exists a unitary Brownian motion (?7")f>o, starting at /, independent 
of Un- Let us define, for each n, 

X" = {Tt[A-,U^UZ], . . . ,TT[AlUnUZ]) 

r« = e-^(Tr[A^f/„],...,Tr[A^[/„]) 

= X"-F" = e-T(TrK[/„(e^f/,'; -/)],..., TrKf/„(e^f/,';-/)]) 

By Step I and Theorem II. 2^ as n tends to infinity, the distribution of converges to the 
one of {Zi tg, . . . , Zk^to)- It follows that for n large enough, 

(28) |E{/(Zi,,„, . . . , Zfc,,,)} - E{/(D")}| < e. 

c) At last, since / is 1-Lipschitz for || ■ ||, for all ra, we have 

|E{/(D")}-E{/(X")}| < E{||D"-X"||} 

< (e-*»E{| Tr[A^[/„]|2 + ■ ■ ■ + | TrKt/JI^})!/^ 

< e-*"/2(^-i Tr[A^v4^*] + ■ ■ ■ + Ti[AlAl*]fl^ , 

the last inequality following from (125]) . By (I27|) . it allows us to claim that for all n, 

(29) |E{/(D")}-E{/(X")}|<e. 
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d) To conclude, note that by the right invariance of the law of f/„, X" and the random 
vector (Tr[A^[/„], . . . ,Tr[A^[/;„]) have the same distribution. By ([26]), ([28]) and ([29]), it 
follows that for n large enough, 

|E{/(Zi, . . . , Zfc)} - E{/(Tr[A-f/„], . . . , Tr[A-f/„])}| < Se. 

It closes the proof of the theorem. □ 
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